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The spin and quadmpole moment of the supermassive black hole at the Galactic center can in 
principle be measured via astrometric monitoring of stars orbiting at milliparsec (mpc) distances, 
allowing tests of general relativistic "no-hair" theorems [f[. One complicating factor is the presence of 
perturbations from other stars, which may induce orbital precession of the same order of magnitude 
as that due to general relativistic effects. The expected number of stars in this region is small 
enough that full iV-body simulations can be carried out. We present the results of a comprehensive 
set of such simulations, which include a post-Newtonian treatment of spin-orbit effects. A number of 
possible models for the distribution of stars and stellar remnants are considered. We find that stellar 
perturbations are likely to obscure the signal due to frame-dragging for stars beyond ~ 0.5 mpc 
from the black hole, while measurement of the quadrupole moment is likely to require observation 
of stars inside ~ 0.2 mpc. A high fraction of stellar remnants, e.g. lOA'f© black holes, in this region 
would make tests of GR problematic at all radii. We discuss the possibility of separating the effects 
of stellar perturbations from those due to GR. 
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I. INTRODUCTION 

The supermassive black hole (SBH) at the center of the 
Milky Way galaxy is surrounded by a compact cluster of 
stars that has been the target of observational surveys 
for more than a decade @, 0, [3, H, fl 0] • Near- infrared 
monitoring of stellar positions using adaptive optics tech- 
niques has allowed orbital reconstruction for roughly 30 
stars at distances ranging from 10° — 10 2 milliparsecs 
(mpc) from the SBH 0, [H, E(| • One of these stars (S2) 
has an orbital period of only ~ 15 yr [HI, [13] and its or- 
bit has been followed for more than one full revolution; 
astrometric data for S2 yield a well-constrained mass for 
the SBH, M. = (3.95±0.06) x 10 6 7W Q (assuming a galac- 
tocentric distance of 8.0 kpc) and a location on the plane 
of the sky that is consistent with that of the radio source 
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The velocity of S2 near periapse is a few percent of the 
speed of light, large enough that relativistic effects like 
advance of the periastron become potentially measurable, 
even on time scales as short as a few years [T3, EM EM H3| ■ 
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No such effects have so far been unambiguously observed 
[l6j |: one complicating factor is the likely presence of a 
distributed mass (stars, stellar remnants, dark matter 
etc.) within S2's orbit which could produce Newtonian 
precession of the same order of magnitude as that due to 
general relativity [U [22j . 

If the SBH is rotating, new phenomena occur for stars 
orbiting at very small separations, r < 1 mpc. Dragging 
of inertial frames and torques from the SBH's quadrupole 
moment Q cause stellar orbital planes to precess, at rates 
that depend respectively on the first and second powers of 
the hole's spin angular momentum J. These spin-related 
effects are small compared with in-plane precession, but 
(in the absence of other non-spherically-symmetric com- 
ponents of the gravitational potential) they contain un- 
ambiguous information about J and Q. In principle, 
observed changes in the orbital orientations of just two 
stars would be sufficient to independently constrain the 
four quantities (J,Q), allowing tests of general relativis- 
tic "no-hair" theorems [l[ . The amplitude of these spin- 
related precessions is very small, of order microarcsec- 
onds (//as) per year as seen from the Earth. Plans are 
being developed to achieve infrared astrometry at this 
level [H,[23]. 

Such measurements will require the presence of at least 
a few bright stars on mpc-scale orbits around the SBH. 
While no such stars have yet been observed, extrapo- 
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lation of the observed stellar densities at distances of 
~ 1 pc from the SBH suggests that of order 10° — 10 2 
stars should be present in this region. Due to their fi- 
nite numbers, these stars will generate a non-spherically- 
symmetric component to the gravitational potential with 
an amplitude that scales as \ ~N m+ , where is the mass 
of a typical star and N is their number. Simple argu- 
ments f CT[j) suggest that such stellar perturbations might 
produce changes in the orbital orientations of test stars 
that are comparable in magnitude to the spin-related ef- 
fects. This would complicate the testing of no-hair the- 
orems by adding what is effectively a source of noise to 
the measured precessions. 

Because the expected number of stars in this region is 
so small, direct iV-body integration of the equations of 
motion for all N stars is feasible. The major technical re- 
quirements are a high degree of accuracy in the iV-body 
integrator and the inclusion of terms describing the rela- 
tivistic accelerations due to the SBH, including spinless, 
spin-orbit, and quadrupole-orbit contributions. 

Here we present the results of a comprehensive set of 
such simulations. Our primary goal is to evaluate the 
degree to which star-star perturbations might obscure the 
signal due to the SBH's spin; hence we focus on changes 
in orbital orientations rather than on the evolution of the 
phase-space variables (r, v) [25|, (26|. We ignore all other 
systematic effects that might limit the ability to carry 
out the high-precision astrometry for stars in crowded 
fiels at the Galactic center [27( • 

In fJTT] we summarize the relevant time scales for or- 
bital evolution near the Milky Way SBH. fjH] presents 
the post-Newtonian ./V-body equations of motion includ- 
ing the lowest-order spin-orbit terms and describes the 
iV-body integrator. Observational and theoretical con- 
straints on the distribution of stars and stellar remnants 
near the Galactic center SBH are summarized in ^IVl 
which also describes the parametrized models used to 
construct the iV-body initial conditions. ifVl summarizes 
the results from the integrations, including estimates of 
the number of stars that can be effectively used to mea- 
sure J and Q. $VI\ discusses how the presence of stellar 
perturbations in the astrometric data can potentially be 
detected and removed from the GR signal. E1VIII sums 
up. 



II. SOURCES OF ORBITAL EVOLUTION 

A. Basic quantities 

The orbital period of a star of semi-major axis a orbit- 
ing around the Milky Way SBH is 



P = , « 1.48a 3/2 yr 



(1) 



where M. is the mass of the SBH and a is the star's semi- 
major axis in units of milli-parsecs (mpc). The second 



relation assumes M. = 4.0 x 10 6 A^o an d m * ^ M,, as- 
sumptions which we adopt in the remainder of the paper. 
The length scale associated with the event horizon of the 
SBH is 



GM. 



1.92 x 10" 4 mpc. 



(2) 



We define x t be the dimensionless spin angular mo- 
mentum vector of the SBH, 



GM. 



0<X<1. 



(3) 



The standard (no-hair) relation between J and the 
quadrupole moment is 



1 J2_ 

c M. ' 



(4) 



We adopt this relation below unless otherwise noted. 

In the regime of interest, stellar orbits around the SBH 
can be approximated as Keplerian ellipses that experi- 
ence gradual changes in their orbital elements, due both 
to the effects of relativity and to perturbations from other 
stars. Here we summarize the relevant sources of evolu- 
tion and their associated time scales under this approxi- 
mation. 



B. Relativistic precession 

1. In-plane precession 

In the orbit-averaged approximation, massless test par- 
ticles orbiting a black hole experience advance of the or- 
bital periapse by an angle [80( 

5vo = A s - 2Ajcosi - ^Ag(l - 3cos 2 i) (5) 

per orbit, where the subscripts S, J, Q denote the effects 
due to the black holes's mass (i.e. the Schwarzschild 
part of the metric), spin and quadrupole moment (the 
Kerr part of the metric) respectively and i is the orbital 
inclination, defined as the angle between the SBH spin 
vector and the stellar orbital angular momentum vector. 
To lowest post-Newtonian PN order, 
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(6a) 



(6b) 



(6c) 
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where e is the orbital eccentricity. Since the 
Schwarzschild contribution exceeds in amplitude the 
spin- and quadrupole contributions to the in-plane pre- 
cession for a(l — e 2 ) > 10~ 4 mpc sw r g , advance of the 
periapse does not contain much useful information about 
the SBH spin Q. 

We define the precession time scale due to the 
Schwarzschild term alone as 



P [yr 



t s 



As(a,e) 
irP{a) 



(7a) 



2 

c a 



a 



e 2 ) 



(7b) 
(7c) 



6 GM. 

» 1.29 x 10 3 yr (l - e 2 ) a 5/2 

The Schwarzschild contribution to the in-plane preces- 
sion is large enough to potentially be detectable via a few 
years' monitoring of identified stars at ~ 10 mpc separa- 
tions from the SBH [jj, [H [H| ■ 

2. Precession of orbital planes 

The gravitational field of a Kerr black hole is not 
spherically symmetric. The dominant non-spherically- 
symmetric effect on test-particle orbits is the coupling 
between the spin of the black hole and the orbital angu- 
lar momentum of the particle, known in the weak-field 
limit as the Lense-Thirring effect (2|| [2t| . Again to low- 
est PN order, the change per orbit of the nodal angle Q, 
isQ 



sn = At- A 



Q COS I 



(8) 



[8l| . Among relativistic effects, precession of orbital 
planes depends only on J and Q. We note that the 
frame dragging contribution in Eq. l[8]). which dominates 
at most distances of interest here, is independent of or- 
bital inclination while the quadrupole-induced precession 
is inclination-dependent. Both precessions leave the in- 
clination with respect to the black hole spin unchanged. 

Defining the associated time scales as in Eq. we 
find 



tj = T- 



P 

4x 



c 2 a(l - e 2 ) 



GM. 



3/2 



1.39 x 10 5 yr (l 



tQ 



P 

3^ 



c 2 a(l 



e 

-i 2 



,3/2 



GM. 

1.34 x 10 7 yr (l 



- 2 ~a 7 ' 2 . 



(9a) 
(9b) 
(9c) 
(9d) 



Figure [T] plots tj and tQ as functions of a and e. 

C. Stellar perturbations 

If there is a star cluster around the SBH, the smooth 
contribution to the gravitational force from the dis- 
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FIG. 1: Time scales associated with precession of orbital 
planes about the Galactic supermassive black hole. tj,tQ\ 
precession time scales due to frame-dragging and to the 
quadrupole torque from a maximally-spinning SBH. Line 
thickness denotes orbital eccentricity, from e = 0.99 (thick- 
est) to e = 0.9 and e = 0.5 (thinnest), tjv: approximate 
precessional time scale due to Newtonian perturbations from 
other stars, assumed to have one Solar mass. Line thickness 
denotes total distributed mass within 1 mpc from the SBH, 
from 10 3 .M© (thickest) to IMq (thinnest), assuming that 
density falls off as r _1 . Shaded (green) region shows range of 
interesting time intervals for observation, 1 yr < At < 10 yr. 



tributed mass breaks the degeneracy between radial and 
angular periods in the classical Kepler problem, causing 
an in-plane precession, in the opposite sense to the rel- 
ativistic periastron advance. Assuming that the stellar 
mass density follows r~ 7 , with r the distance from the 
SBH, the advance of orbital periapse in one period is 



2tt 



M. 



Vl - e 2 F( 7 ) 



(10) 



where M*(r) is the distributed mass enclosed within ra- 
dius r and F = (3/2, 1) for 7 = (0, 1) 0. Setting F « 1, 
the associated time scale is 



^(l_ e 2)-V2 

2 M* 1 ' 



/w = —~rr i 1 -<-) ' - (Ha) 
3.0 x lOVM^a 7 " 372 ^ _ e 2 )- 1 ' 2 (lib) 



where M* is the stellar mass within 1 mpc in units of the 
solar mass. This time scale is long compared with the 
time ts for relativistic periastron advance, Eq. (JT]), at all 
radii of interest unless M* is unphysically large. 
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The discrete nature of the stellar cluster adds an ad- 
ditional, non-spherically-symmetric component to the 
gravitational potential, which can induce precession in 
orbital planes that mimics the effects of frame-dragging 
and quadrupole torques. In the case that the time scale 
associated with this precession is long compared with 
both the radial period (Eq. [1} and with the time scale 
for in- plane periastron advance (Eq.[7|) (assumptions that 
will be verified below), orbits around the SBH respond 
to the finite- TV component of the gravitational force as 
if they were annuli, changing their orientations but not 
their eccentricities ( "vector resonant relaxation" ; [3l[ ) . 

Here we estimate the rate of precession due to finitc- 
N stellar perturbations, adopting a purely Newtonian 
model for star-star interactions. 

Let q = m*/Af # be the ratio between stellar mass and 
SBH mass, N the number of stars and/or stellar rem- 
nants in the region contained within a test star's orbit 
(a more precise definition of N is adopted in §IV) and 
L c the angular momentum of a circular orbit of the same 
energy as that of the test star. In the vector resonant 
relaxation (RR) regime, orbital angular momenta evolve 
approximately as [3l| 



IALI 



(3 v qVN— 



(12) 



for a time At < i co h, where f3 v is a constant of order unity 
and £ co h is the time scale associated with the most rapid 
process that randomizes orbital planes, thus breaking the 
coherence. 

In the absence of GR effects, the only source of 
coherence-breaking is the stellar perturbations them- 
selves ("self-quenching"), for which i co h = tjy, where 
t j\[ is defined by the condition |AL|/L c (tjv) = 1. On 
time scales long compared to fjy, and in the absence 
of frame-dragging or other torques, orbital orientations 
would evolve approximately as 



IALI 



v qVN 



\AcohAt 

p ' 



(13) 



i.e. as (Ai) 1 / 2 rather than as (At) 1 . Fig.[T]shows that t?q 
is > 10 4 yr for reasonable models of the stellar cluster, 
much longer than the ~ 10 yr time scales of interest 
here; hence the self-quenched regime is irrelevant in what 
follows. 

However at some radius, fjy will exceed the time scales 
associated with GR precession of orbital planes, and the 
precession rate will be given by the expressions derived in 
the previous section rather than by Eq. (fT2| . To estimate 
this radius, we begin by expressing the GR precessional 
time scales defined above in terms of the "penetration 
parameter" g = (1 + e)r p /r g > 1, where r p is the Keple- 
rian orbital periapse distance and r g is defined in Eq. 
The results are: 



tQ 



a— 2 P 



;Q X 



(14a) 
(14b) 



1 1 1 1 1 
_m» = M 


I I I I I 
_m, = 1 M 


- ffc^X) - 


- ff^^Y) - 






— 1 — 1 — 1 — 1 — 1 — 


— i — i — i — i — i — 


m, = 3 M Q 


m, = 10 Mq__^^ 


- ( I Sv>^ 1 

\ V \ / J 

1 1 1 1 1 


I I I I I 


-1 -0.5 0.5 1 


-1 -0.5 0.5 1 


L / L o 


L xAo 



FIG. 2: Evolution of orbital planes in a cluster of eight stars 
orbiting about the Galactic center SBH, for an elapsed time of 
2 x 10 6 years. The SBH rotates about the z-axis with maximal 
spin. Four different values were assumed for the stellar masses 
m*, as indicated. Stars were placed initially on orbits with 
semi-major axis 2 mpc and eccentricity 0.5 and with random 
orientations. 



and the vector RR time scale itself is 

1 



tN 



qVN 



P; 



(14c) 



the latter expression is true only up to 0(1) factors which 
have to be derived from simulations. 



Since 



^ = (4/3 X )£ 1/2 >4/3, 
tj 



(15) 



precession due to frame-dragging is everywhere faster 
than precession due to the quadrupole torque. 

The condition that frame dragging dominate stellar 
perturbations is approximately 



In 



F\ 



q 



> i 



(16) 



Fig. [T] shows that for \ = 1, this condition is satisifed 
inside ~ 1 mpc for reasonable values of the enclosed mass. 
For x < 1, the critical radius is smaller. This justifies 
looking at stellar perturbations as a source of "noise" in 
tests of GR. 
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D. Comparing relativistic and Newtonian 
precessions 

Precession of orbital planes induced by stellar pertur- 
bations differs qualitatively from precession due to frame 
dragging since it does not respect the direction of the 
SBH spin axis. Fig.[2]shows the results of a set of 2 x 10 6 
yr integrations (using the algorithm described in ijllip 
that illustrate the difference. When the stellar masses are 
set to zero, orbital angular momenta exhibit the uniform 
precession about the SBH's spin axis associated with 
frame dragging; when stellar masses are increased, the 
orbital angular momentum vectors move quasi-randomly 
about the unit sphere. 

In comparing GR precession with that due to stellar 
perturbations, it is therefore useful to have a measure of 
orientation that is invariant to the direction of the SBH 
spin. We adopt AO, defined as the angle between the 
initial and final orbital angular momentum vectors: 



cos AO = 



L, L 



/ 



LiL 



iff 



(17) 



On time scales of relevance here, |L| is conserved, i.e. 
Li Lf, because the time scales for both non-resonant 
and resonant relaxation are much longer than 10 yr. 

Precession induced by GR changes only the nodal angle 
n (Eq.HJ). Since 

L = L (sinf2sinie x + cosilsinze y + cosie z ) , (18) 

where the z-axis is parallel to the SBH spin vector, Li = 
Lf implies 

cos A9qr — cos 2 i + sin 2 i cos AQ (19) 
which for small Afl is 



A6gr ~ sin iAtt. 



(20) 



In the case of star-star perturbations, we need to ex- 
press A9 in terms of AL/L C . By definition, 



IALI 



L 



2LiL f cosA6> 



so 



|AL| 

LI 



L 2 

2— (l-cosA0) 



(21) 

(22a) 
(22b) 



where the last expression again assumes small AO. 

Specializing Eq. (|2U)) to the case of frame dragging, 
we note that AVL is independent of cosi. Considering 
orbits with a single eccentricity e and with an isotropic 
distribution of inclinations, the rms values of the angles 
in Eq. ([2H)l are therefore related by 



A6j 



An 



(23) 



or (cf. Eq. [6b 



AO j « 4?r 



(24) 



In the case of quadrupole-induced precession, Afl oc 
cosi. Again computing the rms values assuming random 
orientations gives 



A0 Q 



6 2 _ 2 At 

r xg -p- 



(25) 



Finally, for stellar perturbations, we ignore a possi- 
ble dependence of |A0/v| on orbital eccentricity. In an 
isotropic cluster, the orbital angular momenta at any en- 
ergy (~ radius) are distributed as n(L)dL — 2LdL/L 2 , 
so (L 2 ) /L 2 = 1/2, and the rms values in Eq. ([22|) are 
related by 



A9 



N 



V2 



ALI 



Using Eq. (fT2|) . this can be written 

A6 N » V2p v qVN^- . 



(26) 



(27) 



Eilon et al. [32j give f3 v » 1.8 as an average value for 
a cluster with isotropically distributed velocities, if N is 
defined as the number of stars within a sphere of radius 
r = a. We can then write a slightly more accurate defi- 
nition of the vector RR time scale (again defined as the 
time such that A#jv = 



tN 



1 



P 



1.2 



-.P. 



(28) 



This is the expression plotted in Fig. 1. 

We have assumed that precession is due either to 
GR spin effects or to stellar perturbations. In reality, 
one expects vector RR to be quenched somewhat by 
coherence-breaking due to GR precession even at radii 
where ijv < tj,Q- 

Equating (f2"4")l with (|2~7| , we obtain an approximate 
expression for the radius at which frame-dragging domi- 
nates stellar perturbations: 



n 3/2 



4tt x 



V3/3 V q 



r cr it w lmpc(l - e' 



crit 



30 



-1/3 



(io.vi . 



(29) 



-2/3 



.(30) 

where iV cr it is the number of stellar perturbers within 
r cr it of mass m* each. We evaluate this expression in £HVI 
after specifying a model for the stellar distribution, and 
in SjV] we present the results of full iV-body simulations 
that allow more precise estimates of r cr j t . 
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III. A-BODY TREATMENT 

Integrations of the A-body equations of motion were 
carried out using algorithmic regularization [33|, HH im- 
plemented with a chain structure [35| and the time- 
transformed leapfrog [3(| . The algorithm produces exact 
trajectories for Newtonian two-body motion and regular 
results for close encounters involving arbitrary numbers 
of bodies. Velocity-dependent forces were included via a 
generalized mid-point method [37|; the ARCHAIN code [38| 
also incorporates pairwise post-Newtonian forces for non- 

I 



a./,i = 



a J,j = 
X = 

r ij = 



Here, particle number 1 is the SBH and particles j, 2 < 
j < N are the stars. The "linear momentum" that is 
conserved by these equations is 

p = E m ^ + 5 E |5r (x « x Jj) ■ (32) 

i ij % 3 

i 
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IV. MODELS FOR THE STELLAR 
DISTRIBUTION 

A. Observational constraints 

The distribution of stars and stellar remnants at dis- 
tances < 1 pc from the Galactic center SBH is poorly 
understood. Only the brightest stars in the inner parsec 
have been identified, via speckle or adaptive optics imag- 
ing and spectroscopy in the near-infrared bands [ill . |42| | . 
Most of these stars appear to belong to one of two distinct 
populations: (1) "early-type" (ET) stars - apparently 
normal, upper-main-sequence giant stars of O and B 
spectral types with inferred masses of 7 — 80Mq and ages 
less than the main-sequence turnoff age, i.e. OflO 1 — 10 2 ] 
Myr; and (2) "late-type" (LT) stars - old, metal-rich, M, 
K and G-spectral-type giant (post-main-sequence) stars 



spinning particles of orders up to and including PN2.5 
[391 ] . We included PN terms in the interactions between 
the SBH particle and the N - 1 "star" particles. All N 
particles were included at all times in the chain. Accu- 
mulated energy errors were never more than a few parts 
in 10 10 . 

We modified ARCHAIN to include the lowest-order con- 
tributions of the SBH's spin and quadrupole moment to 
the motions of the stars. In the covariant spin supple- 
mentary condition (SSC) gauge [4(|, the spin-related, TV- 
body accelerations aj are 



(31a) 

(31b) 

(31c) 
(31d) 



Adopting Eq. ^ for the SBH quadrupole moment, the 
equation of motion for the jth particle has the additional 
term slqj, where 



J=J/J. (33) 

I 

with ages O[10° - 10 1 ] Gyr and masses 1 - 2M Q . The 
density of ET stars increases steeply toward the SBH and 
these stars account for a large part of the total luminosity 
of the central cluster, but their total numbers are small, 
roughly 10 2 in the inner 0.1 pc [(| [13, IHI with few if 
any on orbits that bring them within ~ 10 mpc from 
the SBH, making them unlikely candidates either as test 
stars for observing GR spin effects or as perturbers of the 
test stars. 

The LT stars on the other hand are believed to be 
characteristic of the dominant, old population; roughly 
6000 LT stars have been identified in the inner ~ 0.5 pc 
and their K-band luminosity function suggests a roughly 
continuous star formation history over the last ~ 10 Gyr 
[H^l^. In spite of their large numbers, the observed LT 
stars appear to be weakly concentrated toward the SBH. 
Number counts complete to K sa 15.5 (corresponding 
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to the sub-giant phase for 1A4© stars) reveal a projected 
density that is flat or declining insid e a p rojected distance 
of ~ 0.5 pc from the SBH 03, M, El- while the exis- 
tence of four LT stars on very tight (5 mpc < a < 20 mpc) 
orbits around the SBH has been established [10| , depro- 
jcction of the binned surface density profile implies a cen- 
tral space density that is consistent with zero at distances 
smaller than ~ 0.1 pc from the SBH [471 ] . 

The low density of LT stars in the inner parsec is not 
well understood. If the time scale for exchange of orbital 
kinetic energy between stars (the two-body relaxation 
time; 48]) is shorter than several Gyr, one expects the 
stellar distribution to have attained a quasi-steady-state 
distribution of the form n(r) ~ r r , 3/2 < 7 < 7/4 
[49I [50( ] within the SBH gravitational influence radius, 
r, n fl = GM m /a1 ~ 10° pc. This is clearly not observed 
[45j . suggesting either that the relaxation time exceeds 
~ 10 Gyr throughout the inner parsec, or that the bright- 
est stars have been hidden from view or destroyed. Col- 
lisions with main-sequence stars or stellar remnants can 
remove the outer envelopes of red-giant stars, pot entially 
explaining the low observed density of giants [42| . How- 
ever this mechanism only appears to be effective at dis- 
tances less than ~ 0.1 pc from the SBH [HT|, [52|, [H| , even 
assuming a high density for the colliding populations (an 
assumption for which there is currently no observational 
support). Even at these radii, collisions would seem to 
be ineffective at explaining the depletion of stars down 
to magnitudes as faint as 15.5 [53| . 

It has been argued that the stellar initial mass function 
(IMF) may have been strongly truncated below ~ 3Mq 
in the Galactic center region [54|, [55[ . These are just the 
stars that would dominate the K-band number counts 
now [5J]. 

At the high-mass end, standard IMFs [H, H3] predict 
that ~ 0.1% of stars have initial masses greater than 
20.M©, ending their short lives as ~ 5 — 15.M© black 
holes (BHs). The BHs are expected to segregate nearer 
to the SBH than the lower-mass components (stars, white 
dwarves, neutron stars [HI), possibly domina ting the 
total number density inside ~ 1 mpc [52l l59l . l6dl ] and 
providing the bulk of the perturbations acting on the 
observed stars in this region. However if the observed 
distribution of late-type stars is a guide, the two-body 
relaxation time may be too l ong for establishment of a 
mass-segregated distribution [U HH, [6l| . If this is the 
case, there is no compelling reason to assume that the 
ratio of BHs to stars is as large as implied by the mass- 
segregated models. 

Proper motion studies of large samples of LT stars in 
the inner parsec [62l . [63| yield dynamical constraints on 
the distributed mass (stars, stellar remnants, gas etc.) in 
this region. The proper motion data robustly require an 
extended mass of ~ 0.5— 1.5 x 10 6 A4© within the central 
parsec [63] . However these data do not strongly constrain 
the radial dependence of the distributed mass density nor 
the amount of mass on the mpc scales of interest here. 



TABLE I: Parameters of the TV-body models 
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1 
1 
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2 
2 
2 
2 
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B. Parametrized models 

Given these uncertainties, we explored a range of dif- 
ferent models for the distribution of stars and stellar rem- 
nants near the Galactic SBH. We define M+ as the dis- 
tributed mass within 1 mpc from the SBH and M+ = 
M^/Ai®. We idealize the stellar populations in this re- 
gion as consisting of just two components: 1.M© main- 
sequence (MS) stars and 10.M© BHs. The first popu- 
lation is assumed to be amenable to astrometric moni- 
toring, and all discussions of orbital evolution presented 
below will refer to this population. While the orbits of 
the BHs are also allowed to evolve in our models, we do 
not describe that evolution in what follows. 

In addition to M + , three additional parameters define 
the initial distributions of stars and stellar remnants in 
our models: 

• the power-law index 7 describing the number den- 
sity profiles, n(r) cx r 7 is assumed to be the 
same for both MS stars and BHs; 

• the (number) ratio 1Z of BHs to MS stars, i.e. 1Z — 
Nbh/N ms ; 

• the velocity anisotropy /?, defined such that 
of /of = (1 — 1 where a r and a t are respec- 
tively the Id velocity dispersions in directions par- 
allel and perpendicular to the radius vector; (3 = 
corresponds to isotropy. 

For M+ we adopted one of the three values 
(10, 30, 100)A4©; the latter value is roughly the enclosed 
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FIG. 3: Approximate value of the radius at which stellar per- 
turbations match frame dragging in terms of their ability to 
change the direction of orbital angular momenta (Eq. [29}. 
Two models for the stellar cluster are shown: a steeply-rising 
density profile, 7 = 2, that is dominated by IOMq stellar 
black holes; and a shallower density profile, 7 = 1, dominated 
by 1A4q main-sequence stars. Horizontal axis is the total 
distributed mass within 1 mpc from the SBH. Line widths 
denote SBH spin: \ = 1 (thickest), \ = 0.3, and \ = 0.1 
(thinnest). 



mass predicted by the relaxed, mass-segregated models 
cited above. For 7 we considered the values (0,1,2); 
7 = corresponds to a constant density in the inner 
mpc, roughly what is observed in the projected density 
of LT stars, 82] while 7 — 2 is approximately the value ex- 
pected for a mass-segregated population around a SBH. 
Adopting the proper-motion result that the distributed 
mass within 1 pc is ~ 10 6 A4q [H|, the implied mass in- 
side 1 mpc is ~ 10°(10 3 ).M Q for 7 = (1, 2). We therefore 
associated larger values of M* with larger values of 7, 
although as noted above, the proper-motion data do not 
directly constrain the mass distribution on mpc scales. 
TZ was set to (0,0.1,1); TZ — 1 is roughly the largest 
value predicted in the mass-segregated models assuming 
a standard IMF, while TZ m 10 -3 is expected in the ab- 
sence of any mass segregation. For the small (~ 10 2 ) 
total particle numbers in the A^-body simulations, 7Z = 
is essentially the same as TZ = 10 -3 . 

A steady-state orbital distribution in a point-mass po- 
tential requires f3 < 7— 1/2, i.e., isotropic velocity distri- 
butions are not permitted when 7 < 0.5: the distribution 
of orbital eccentricities must be biased toward small val- 
ues when the spatial distribution is flat. The Galactic 
center proper motion data cited above [H[ suggest ap- 
proximate isotropy in the (projected) inner parsec. The- 
oretically, two-body encounters should drive the distri- 
bution toward isotropy near the SMBh while at the same 
time populating the low-angular-momentum orbits, pro- 



ducing an isotropic density cusp. However since the cusp 
is not observed, it is not clear that relaxation has had 
sufficient time to reduce anistropies to low values [47j . 
We therefore considered non-zero values of even when 
setting 7 = 1 or 2. 

The following distribution of orbital elements: 



N(a, e 2 )dade 2 

g( a ) 

He 2 ) 



N g(a)h(e 2 )dade 2 



- „2- 7 



(34a) 
(34b) 



1 



e 2 ) , (3 < 7 - 1/2 (34c) 



generates steady-state phase-space distributions with the 
properties defined above. Monte-Carlo realizations of the 
stellar positions and velocities were generated from this 
expression given the parameters (7, /?, M*, TZ). Relativis- 
tic corrections were ignored when generating the initial 
conditions. We assumed that N(a, e 2 ) = for a > a max ; 
in most of the simulations, a max = 4 mpc, but smaller 
values were adopted as necessary to limit the total num- 
ber of particles to ~ 180, since for larger N the ARCHAIN 
routine was found to run very slowly. Table U gives the 
composition of all the models discussed below. 

Stellar orbits were excluded from the initial conditions 
if their periapse fell below 20GM./c 2 , roughly the dis- 
tance at which a solar-mass star on the main sequence 
would be tidally disrupted. 

For each set of parameters defining the initial distri- 
bution, a set of different Monte-Carlo realizations was 
generated and independently integrated forward in time. 
The number N iaa n of independent realizations was chosen 
such that the total number of MS stars in the combined 
set of integrations was > 10 3 . In addition, some inte- 
grations were repeated with the masses of the MS stars 
and BHs set to zero (leaving just the relativistic terms 
capable of inducing evolution of the orbital planes); and 
with the SBH spin set to zero (leaving just the Newtonian 
perturbations). 

Given these models for the stellar cluster, we can use 
Eq. (|29p to compute the approximate radii r cr ;t where 
frame-dragging begins to dominate stellar perturbations. 
Assuming that the density is dominated either by IMq 
MS stars (TZ < 1) or by lftM© BHs (K > 1), we find: 



^*crit 



(43, 14) X 



(1 - e 2 ) 3 / 2 M* 



1/2 



mpc 



(35) 



where the first number in parentheses refers to the MS 
cluster and the second to the BH cluster. Fig. [3] plots 
r cr i t vs. M+ for various values of x assuming e = 2/3, 
the mean eccentricity in an isotropic distribution. 



V. RESULTS 

In the A^-body simulations, we characterized changes 
in stellar orbital orientations in two ways: via Afi, 
the change in the nodal angle (defined with respect 
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FIG. 4: Changes over 10 years in the orientations of stellar orbital planes, as measured via AS1 (top) and A9 (bottom). 
Parameters of the iV-body models were 7 = 2, /3 = 0, H = 1, and M* = 30.M© (Nms = 10, JVbh = 11). Each point 
corresponds to a single star in a single integration; red points are orbits with initial eccentricities 0.7 < e < 1 and blue points 
have < e < 0.7. In the lower panels, dashed lines show Eq. (I24|l . the frame-dragging precession, for e = 2/3, and dotted lines 
show Eq. (|27[) . the approximate model for precession due to stellar perturbations. 



to the SBH equatorial plane); and via the coordinate- 
independent quantity A0 7 the angle between the initial 
and final orbital angular momentum vectors (Eq. I17[) . 
The nodal angle advances uniformly in time in response 
to GR effects (Eq. [5]); furthermore for a given \ an d Ai, 

the quantity (l — e 2 ) 3 ^ 2 AJ7 depends only on a in the 
frame-dragging regime (Eq. [6bl [8| • 

Fig. H plots (1 - e 2 ) 3/2 An and AO vs. a for each of 
the MS stars in a set of 10-year integrations of models 
with 7 = 2, (3 = 0, TZ = 1 and M* = 30 and three differ- 
ent values of the SBH spin, x = (1,0.1,0). Also shown 
are the predicted, rms values of AO from Eqs. (f2~4")) (frame 
dragging) and (|27|) (stellar perturbations). In this model 
cluster, in which stellar BHs dominate the total mass, 
stellar perturbations dominate changes due to frame- 
dragging beyond radii of ~ 1(0.3) mpc for x = 1(0.1). 
As measured via AO, GR effects are strongest for eccen- 
tric orbits, as expected, while the amplitude of the stellar 
perturbations is not noticeably e-dependent. The stellar 
perturbation model derived above is reasonably good at 
predicting the mean value of AO in the integration with 
X = 0, although the observed dependence on a appears 
to be shallower than predicted for a < 0.5 mpc. 

Because of the large scatter in the amplitude of stel- 
lar perturbations at each a, the radius at which the GR 
signal clearly stands out from the "noise" is somewhat 
smaller than would be predicted from the rms values 
alone (Fig. [2]). 

The dependence of these results on the amount of 
distributed mass is shown in Fig. which summa- 



rizes results from integrations of models with M* 
Mil. 30, lOO)^©; other parameters are as in Fig. 0] As 
is increased, the amplitude of the "noise" from star- 
star perturbations increases, roughly in proportion to 
Af*. These plots also indicate the expected amplitude 
of the quadrupole-induced precession. For x = 1, stellar 
perturbations dominate changes due to the quadrupole 
at radii beyond ~ 0.5(0.3) mpc for M* = 10(100)A^ Q . 

Fig. [5] shows the results of a comprehensive set of 
integrations using different models for the stellar clus- 
ter (Table 1). Especially when the SBH spin is low 
(x = 0.1), stellar perturbations can dominate the signal 
due to frame-dragging down to very small distances from 
the SBH, e.g. ~ 0.2 mpc for M+ = 100, corresponding to 
orbital periods of ~ 0.1 yr. 

We define the "astrometric precessions" AG = 
(a/D)A0 where D — 8.0 kpc is the distance to the Galac- 
tic center 

Q- AG is roughly the angular displacement of 
the orbital axes, as seen from the Earth (ignoring projec- 
tion effects). 

In the following section we discuss how measurements 
of J and Q might be feasible even in cases where the 
stellar perturbations are significant. Here, we assume 
that, in order to be useful for tests of GR, a star must 
satisfy two minimum conditions: 

1 Its astrometric precession must exceed some mini- 
mum threshold set by the detector. 

2 Its precession must be dominated by GR effects. 

We call MS stars that satisfy both conditions "de- 
tectable." 
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FIG. 5: Similar to Fig. [4] except that average values have been computed in bins of semi-major axis. Three different JV-body 
models are shown, differing in the distributed mass: M* = (10, 30, 100)^©. All models have 7 = 2, /3 = 0, 1Z — 1 as in Fig. [4] 
Open circles: \ = 1; squares: x = 0.1; crosses: \ = 0. The predicted angular changes for \ = 1 due to frame-dragging are 
shown as the dashed lines in the upper histograms and as the yellow band in the lower histograms. The cross-hatched regions 
indicate the range of precession amplitudes expected from the SBH quadrupole moment alone. Dotted lines in the lower frames 
are the approximate model for stellar perturbations, Eq. (|27l) . 



We base our assumptions about the minimum observ- 
able angular changes on the specifications for the planned 
instrument GRAVITY GRAVITY will observe the 
Galactic center three times a year (in April, July and 
September) and the error in each astrometric data point 
will be 5Q ~ 10/ias (F. Eisenhauer, private communica- 
tion). Since the precession is linear with respect to time, 
the uncertainty in the measured astrometric precession 
after n observations is 

cr A e~2v / 3W " 1 56«35/iasW ™ 1 . (36) 
y n{n + 1) y n(n + 1) 

For an elapsed time of 1(3)10 yr, ctao = 14(10)6. 2/ias. 

Fig.[7]plots the distribution of A6 values after 10 years 
due to star-star perturbations in one model, (7, /?, M*) = 
(2,0,30). Also plotted are the distributions that would 
arise from frame-dragging and quadrupole torques alone, 
Ae^g = {a/D)Mlj^ Q sini (cf. Eq. (201), for the same 
stars. Fig. [7] suggests that a clean separation of frame- 
dragging and Newtonian precessions for most stars in 
this model requires a < 0.5 mpc for x — 1 an d a < 0.2 
mpc for x — 0.1. For these values of a the amplitude of 
A8j is greater than 10/xas for most stars, making them 
accessible to astrometric monitoring. High eccentricities, 
e > 0.9, allow these requirements to be relaxed: eccentric 
orbits with a as large as ~ 1(0.5) mpc can produce mea- 
surable displacements that are signficantly greater than 
those due to stellar perturbations. 

Detecting the effects of the quadrupole moment in this 
model would be considerably harder. For x — lj the 



quadrupole precessions separate cleanly from the stellar 
perturbations only for a < 0.1 mpc, or a < 0.3 for the 
most eccentric stars. 

Regardless of the number of stars that satisfy condi- 
tions 1 and 2, an additional requirement is that 

3 the stars satisfying these conditions must consti- 
tute a large fraction of all stars in the region being 
observed 

since otherwise there is a large probability that the pre- 
cession of a randomly-chosen star will be dominated by 
non-GR effects. 

We calculated the average number of detectable stars, 
(TVdetect), in each of our models after expressing condi- 
tions 1 and 2 in the forms 

1 A6 > cr A e 

2 A6 JtQ > A6 95 

where A9gs is the upper edge of the 95% confidence in- 
terval of the distribution of stellar perturbations. Obser- 
vational intervals of At — (1, 3, 10) yr were considered. 

The results are shown in Figs. [5] and [5] for four models 
of the stellar cluster, and in four bins of semi-major axis. 
The figures also show /, the ratio of (Adotcct) to the total 
number of stars in each bin, and ( Af2) , the average value 
of the "astrometric precession" for the detectable stars. 
(iVdetect) ^ 1 means that at least one star would be ex- 
pected to be present that satisfies the two detectability 
criteria; if in addition / is large, such stars constitute a 
large fraction of all stars in the same radial bin. 
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FIG. 6: Changes over 10 years in the orientation of stellar orbital angular momentum vectors, in JV-body integrations of various 
models (Table 1). The three columns correspond to three values 1Z = (0, 0.1, 1) of the ratio of BHs to MS stars; the values of 
7, f3 and M* that characterize the stellar distribution are given in the lower left of each panel. Circles are for integrations with 
X = 1 and squares are for \ = 0.1. Dashed lines show the expected, rms contribution to A8 from frame dragging for both 
values of %■ 
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The figures illustrate the tradeoff that occurs between 
the number of detectable stars at a given radius, and the 
certainty that a single star at that radius is responding 
to GR effects rather than to stellar perturbations. For 
example, going across the second column in Fig. [8] (i.e. 
7 = 1,0= 0,M* = 30), (N) increases steadily with 
increasing a, but / behaves oppositely, dropping below 
10% in most cases for a > 1 mpc. 

With regard to frame dragging Fig. [8] suggests the fol- 
lowing. 

• In models with low central densities, 7 = (0, 1), de- 
tection of frame-dragging precession may be feasi- 
ble after At > 3 yr for orbits with 0.2 < a/mpc < 1. 
At smaller radii the number of stars is too small; 
at larger radii the noise from stellar perturbations 
is too great. 

• In models with a steep density profile, 7 = 2, detec- 
tion of frame dragging is feasible at a < 0.2 mpc af- 
ter At — 1 yr in most of the models; at a < 0.5 mpc 
after At = 3 yr; and at a < 1 mpc after At = 10 
yr. The exceptions are models with a large fraction 
of stellar BHs (TZ—1) in which the stellar perturba- 
tions always dominate. 

With regard to quadrupole precession, Fig. [9] sug- 
gests that (-/V<jetect) > 1 occurs in tandem with large 
(> 50%) / only for rather narrow sets of parameters, e.g. 
7 « 2,M « 30, a < 0.2mpc,ft < 0.1, At > 10 yr. De- 
tecting the effects of the quadrupole torque above noise 
from the stellar perturbations is apparently only feasible 
if the stellar cluster is rather finely tuned: there must 
be a substantial number of stars very close to the SBH, 
r < 0.2 mpc, but at the same time a small number of 
stellar remnants so that the stellar perturbations do not 
dominate. 

In practice, the quadrupole-induced precession may be 
large compared with that due to stellar perturbations, 
but still small compared with precession due to frame 
dragging, making it difficult to see in the data. As shown 
in Fig.[TJ £q < tj only holds at very small radii, r < 0.05 
mpc, unless eccentricities are large. 

VI. DISCUSSION 
A. Compensating for stellar perturbations 

The numerical experiments described above were de- 
signed to elucidate the extent to which "noise" from 
stellar perturbations can mask the signal from GR spin 
effects at the Galactic center. In situations where the 
stellar perturbations are present but not dominant, one 
would like to be able to detect the perturbations and re- 
move their effects from the data. Here we outline one 
approach to that problem. 

Suppose that one observes a set of stars at the discrete 
times tj,j — l,...,7V bs- Denoting the stars by index 
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FIG. 8: Detectability of frame-dragging in four models for the stellar cluster. Results are displayed in four bins of semi-major 
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angle of the detectable stars in the bin, expressed in units of 10~ 4 arc min. In each group of similar symbols, the left (red), 
middle (blue) and right (black) symbols are for models with 1Z = (1, 0.1, 0) respectively, while the integration time is indicated 
by the size of the symbol: 1 yr (smallest), 3 yr, and 10 yr (biggest). The tick marks indicate (-/Vdctcct) = 1, / = 50%, and 
(AG) = 10/xas. 
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k = 1, . . . , iV or b, the A^rb orbital solutions are required 
to satisfy the equations 

= u k x L k (tj) + r(tj) x F N (tj,r) , (37a) 
u k = P^iAj^- A Qtk cos i k ) J (37b) 
Here Fjv(i, r) is the perturbing force field due to all the 



stars and stellar remnants, and r k {t) is the position of the 
fcth star at time tj. Fjy depends explicitly on the time, 
to the extent that the orbits of the stars that produce 
the perturbing force are changing. 

We first consider the orbital solutions for each ob- 
served star separately. In their most general form, equa- 
tions (|37[) are under-determined (3-/V b s equations with 
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3 + 3A^ bs unknowns). To make progress, it is necessary 
to introduce approximations. On the short time scales 
of relevance, the contribution to the total gravitational 
potential from the stars can be considered (1) fixed in 
time and (2) not strongly varying in space. For instance, 
a spatially uniform F is the lowest order term in a mul- 
tipole expansion, and one could add, if necessary, dipole 
and quadrupole terms, etc. 

As an example, we consider a stellar perturbation that 
is modeled as 



$*(x, y, z) = 2irG Pt (A x x 2 + A y y 2 + A z 



(38) 



This is the potential due to a homogeneous triaxial clus- 
ter centered on the SBH and aligned with the coordinate 
axes; the dimensionless quantities (A x , A„, A z ) are deter- 
mined by the axis ratios of the ellipsoid [6J] . Proceeding 
as in the derivation of Eq. ((SJ , one can derive the (New- 
tonian) orbit-averaged effect of this perturbing potential 
on the orientations of orbital planes. For instance, in the 
case of low-eccentricity orbits, one finds 65 1 



dL 

dt 



-N 



4nGp t a 



2 ( {A z -A x )L y L z 
{Ay — A Z )L X L Z 



(A x 



(39a) 
(39b) 



Ay)L x Ly 



with L constant. Similar expressions can be derived for 
orbits of arbitrary eccentricity [651 ] . 

In practice, one would simultaneously fit all the pa- 
rameters appearing in Eq. (|3T[) to the observed positions 
and/or velocities. In the case of the model potential just 
assumed, this would amount to introducing 9 additional 
parameters to the 6 required for a Keplerian solution, i.e. 
ui)-; the perturbing density p t and its two axis ratios; and 
the three direction cosines that define the orientation of 
the ellipsoid. A data set with N ^ s > 5 is then formally 
well-determined. Note however that iV bs ^ 5 is required 
to obtain a full solution with a reasonable goodness-of- 
fit. Repeating this fitting procedure for the N or b orbits 
will then allow one to estimate the orbit-to-orbit varia- 
tion in the parameters that define the perturbation. If 
the latter is large, the model can be made more general 
by the addition of extra parameters. 

This fitting procedure will yield an estimate of uik for 
the iV or b orbits that is independent of any assumptions 
about the underlying physics of the GR precession, apart 
from the assumption that the precession rate is constant. 
Given enough orbits, the empirically determined values 
ujk can be correlated with the orbital properties to test in 
an almost assumption-free way non-standard theories of 
gravity. The no-hair conjecture can be tested explicitly 
be substituting into equation (|37b[) the GR expressions 

0: 



3/2 



Aj.k = 47TX0J. 
A.Q,k = 3nx 2 Qk 2 



with q the penetration parameter defined above. The 
torque equation (|37[) can then be written as 



dLfe / -3/2 ,,> 

l + ^p- 1/2 X-(U/L k ) 



(xxL fc )+T„, fc (41) 



The second (small) term in the square brackets (pk ^S> 1) 
reflects the relative contribution of the quadrupole mo- 
ment to the precession. The same fitting procedure out- 
lined above for the general case can be carried out here, 
with the difference that the vector %, which replaces Wk, 
is common to all the orbits, and therefore a simultane- 
ous fitting of the entire set of iV or b orbits will improve 
the power of the orbital solutions. The comparison of 
the best fit values of x from individual orbits to the one 
from the simultaneous global fit can help assess the ro- 
bustness of the result. 

An initial test of the no-hair conjecture could be ob- 
tained by repeating the fit once with, and once without 
the quadrupole term, to check whether the quadrupole 
term indeed improves the fit. A slightly more discrimi- 
nating test would be to introduce an extra free pre-factor 
Jq to the quadrupole term, and see whether the best 
fit yields /q ~ 1. More sophisticated tests would re- 
quire assuming specific alternative functional forms for 
the quadrupole term. 



B. Independent constraints on the spin 

The foregoing assumes that the mass M. of the Milky 
Way SBH is a known quantity, i.e. that it need not be 
treated as a free parameter when fitting the astrometric 
data to Eq. (|37[) . Determinations of M, are based on 
the motion of stars on spatial scales of 10 1 — 10 3 mpc 
SOS [HI, outside the region where GR spin effects are 
detectable. 

In the same way, independent measurements of the 
magnitude and direction of the SBH's spin J could 
remove as many as three additional parameters from 
Eq. (|37[) . reducing the uncertainty on estimates of the 
quadrupole moment. 

A number of approaches to the determination of J are 
being pursued. Sgr A*, the compact source of radio, in- 
frared and X-radiation at the center of the Milky Way, 
exhibits variability at IR and X-ray wavelengths on time 
scales as short as 20 — 30 minutes that can be interpreted 
as emission from "hot spots" in orbits just outside the 
SBH event horizon [6(1 l67l . l68l . [69| . If this interpreta- 
tion is correct, x <^ 0.3 — 0.5 is required in order that 
the orbital period at the innermost stable orbit be as 
short as the observed rise times for the flaring events 
[IS IzOj IzH- Very long baseline interferometry at sub- 
mm wavelengths has the capacity to resolve the structure 
of the plasma surrounding the SBH on angular scales of 



(40) 10 — 100 //as, comparable to r g [72j,l23j. Such observa- 



16 



tions can constrain the spin via comparison with theoret- 
ical accretion disk models, or by detection of variability 
on spatial scales of ~ r g 0, [la, [lH| . The properties of 
the X-ray polarization from the inner accretion disk are 
also predicted to be spin-dependent fri\ . 

Constraints on J derived from observations like these 
will be highly model-dependent, requiring assumptions 
about the nature of the emission, the geometry and phys- 
ical state of the emitting gas (accretion vs. outflow or 
jets), etc. It is not clear at present whether the sys- 
tematic uncertainties of these independent spin measure- 
ments will be larger or smaller than the uncertainties as- 
sociated with astrometric estimates of the spin; the latter 
are due both to the stellar perturbations modelled here, 
and also to the difficulties associated with measurement 
of proper motions of faint stars in crowded fields [13] ■ Ad- 
ditional, if uncertain, information on the SBH spin can be 
used to improve the global orbital solutions by including 
it as prior probabilities in Bayesian best-fit methods for 
large-dimensional parameter spaces, such as the Markov 
Chain Monte Carlo method [Za, [7§| . 

VII. CONCLUSIONS 

1. The spin and quadrupole moment of the Galactic 
center supermassive black hole (SBH) can in principle 
be measured by observing the precession of the orbital 
planes of stars in the inner milliparsec (mpc) [l| . How- 
ever, gravitational interactions between stars in this re- 
gion are likely to induce orbital precession of the same 
approximate amplitude as the precession due to frame 
dragging. The stellar perturbations manifest themselves 
as coherent torques over short time scales, mimicking 
general relativistic (GR) precession. 

2. The number of stars and stellar remnants (e.g. stellar- 
mass black holes, BHs) in this region is uncertain, but 
small enough (~ 10° — 10 3 ) that full Y-body simulations 
are feasible. A regularized post-Newtonian V-body al- 
gorithm is presented that includes the lowest-order spin- 



orbit and quadrupole-orbit terms. 

3. Assuming near-maximal spin for the Milky Way SBH, 
detection of frame-dragging precession may be feasible 
after a few years' monitoring with an instrument like 
GRAVITY (23| for orbits in the radial range between 
~ 0.2 mpc and ~ 1 mpc. At smaller radii the number 
of stars is too small, while at larger radii the star-star 
and star-remnant perturbations dominate GR effects. In 
models where the number of stellar BHs is comparable 
to the number of observable stars, GR effects are almost 
always swamped by perturbations from the remnants. 

4. Quadrupole-induced precession stands out clearly 
from stellar perturbations only in a narrow class of mod- 
els for the nuclear star cluster, having moderate to high 
central densities and a small BH fraction, and only at 
radii r < 0.2 mpc. 

5. Because the orbit-averaged torques from stars are ap- 
proximately constant in magnitude over year-long time 
scales, it is possible in principle to disentangle the effects 
of stellar perturbations from those due to GR, allowing 
tests of gravity even in the presence of stellar perturba- 
tions. 
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